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Abstract. Wolf-Rayet (WR) stars near solar metallicity are believed to be driven by radiation pressure on the UV spectral 
lines of metal ions. As the metallicity decreases so does the line opacity, therefore the mass-loss rate. However, since the 
composition of a WR atmosphere is determined by the burn products of the core, there is a lower limit on the line opacity - 
and therefore the mass-loss rate - of a WR star, even in a star with zero initial metallicity. This presentation is the result of 
attempt to calculate the mass-loss rate of a Population Ill-type WO star using a modified version of the CAK approximation. I 
find that n e > 10 13 cm" 3 and 0.5 < T < 0.7 give the most plausible results, with the resulting mass-loss rate between 2 x 10~ 9 
M yr~' and 3 x 10~ 8 M yr _1 . 
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INTRODUCTION 

At low metallicity line-driven Wolf-Rayet (WR) winds 
exhibit smaller mass-loss rates than WR stars near solar 
metallicity. Since WR winds are significant contributers 
of chemical enrichment and energy deposition at all 
metallicities, understanding the physics of wind driving 
at low metallicity is important. 

WR winds have been studied computationally via ra- 
diative transfer [1] or radiation hydrodynamics simula- 
tions 0. As a complement to these analyses, this study 
uses simple analytic models based on a modified version 
of CAK theory [i3| to determine the effects of a Popula- 
tion III abundance profile on WR wind driving. To keep 
the results as general as possible, model-specific param- 
eters such as the velocity law and the location of the crit- 
ical point are not specified. 

Recently Vink and de Koter H explored the metallic- 
ity dependence of WR mass loss using a Monte Carlo 
approach with CAK-type line driving for stars at metal- 
licity down to Z/Zq = 10~ 5 . They found that the mass- 
loss rate flattens at low metallicity (Z/Z < 10~ 3 ) be- 
cause the self-enrichment of N and He in WN stars and 
C in WC stars provide the lines needed for wind driv- 
ing. They found M « 10~ 8 M Q yr~ l for low metallicity 
WN stars and M = 1.4 x 10~ 7 M yr _1 for low metal- 
licity WC stars. With these results in mind, this presenta- 
tion provides results for a WO star at the extreme of zero 
background metallicity, varying the luminosity (via the 
Eddingtom parameter F) and the electron density. 

CAK THEORY AND MODIFICATIONS 

The first successful analytic models of line-driven winds 
were developed by Castor et al. [3| (hereafter CAK) 
to model the winds of OB stars. CAK theory models 



wind driving via radiation pressure on a large number 
of UV lines. It relies on the Sobolev approximation Q, 
which assumes a rapidly accelerating wind such that a 
photon's Sobolev length L so ^ — v,^(dv / 'dr)~ l , where v,/, 
is the thermal speed and dv/dr is the gradient of the 
wind velocity, is much smaller than the mean free path 
between lines. 

In the CAK model the radiation pressure on the lines 
is expressed as a ratio of the radiative acceleration of the 
lines gi to the radiative acceleration of the free electrons 
g e using the so-called force multiplier M(t), 

M(t) = ^. (1) 

ge 

The force multiplier is then parameterized by the relation 

M(t) = kr a , (2) 

where k and a are determined from the line strength 
distribution, and a in particular is a number between 
and 1 that is related to the fraction of optically thick lines, 
thus a = 1 means all lines are optically thick over L so t,. 

The details of the model used in this analysis, which 
include the effects of frequency redistribution due to line 
branching and bound-free opacity, can be found in Onifer 
and Gayley |6| (hereafter OG) and references therein. 
The following are some highlights. 

In CAK theory the momentum equation at the critical 
point in units of the effective gravity can be written (see 
OG, Sec. 2 and Eq. 27), 

1 +y C = Y~f M (') F NW Fred, O) 

where the first term on the left-hand side represents the 
effective gravity, Y is the Eddington parameter, F^w is a 
correction for multiple scattering (7), F ret j is a correction 
for frequency redistribution due to line branching and 



-free opacity, and y is the inertia scaled to the 
ve pravitv. 

c v c (dv/dr) rc 



bound- 
effective gravity 



y c = 



GAf.(i-r) ' 



where r c and v c are the radius and velocity at the critical 
point and M* is the stellar mass. 

MODEL SETUP 

One cannot hope to model analytically all the complexity 
that a radiative tranfer or radiation hydrodynamics code 
is capable of modeling, so some simplifying assumptions 
must be made. The wind is assumed to be spherically 
symmetric, smooth on scales larger than L so b, and time 
steady, and the star is assumed to be nonrotating. 

At low Z a WR star is more likely to evolve to the WO 
phase than at high Z because the small amount of mass 
lost in the WN phase allows more production of O via 
a-particle capture of C before the overlying envelope is 
stripped. The WO phase can be defined by the relation 
(Nc+No) I > 1-0, where Nj is surface number frac- 
tion 0. The surface abundances are based on the SMC 
WO star Sand 1 studied by Kingsburgh et al. J51. The 
mass fractions of the included elements are Y = 0.21, 
X c = 0.51, X = 0.25, X Ne = 0.02. The Ne is entirely 
22 Ne, which is a product of the He-burning process. All 
ion states of each element are included. The Sobolev 
opacities are calculated using oscillator strengths from 
the Opacity Project [10]. Details of the opacity calcula- 
tion can be found in OG. 

The mass of the star is 5 M , also based on Sand 1. 
This is may be a bit small even for a highly-evolved Pop- 
ulation III WR star due to the lower mass-loss rates of 
such stars, so the mass-loss rates should be thought of 
as lower limits. The terminal speed Vj n f is set such that 
vm — v c = 4700 km s . Kingsburgh et al. [9| report 
Vj n f = 4200 km s _1 , which gives v c = 500 km s _1 . This 
reflects the large terminal speeds of WO stars. M, though 
ultimately has a small dependence on v; n f, since M is set 
by conditions near the surface at the critical point and 
Vj n f is set by the global wind conditions. The tempera- 
ture at the critical point T c = 1.3 x 10 5 K. Two sets of 
models are run. One set assumes complete frequency re- 
distribution (CRD). The other relaxes the CRD assump- 
tion by accounting for for changes in the wind ionization 
structure due to the changing wind temperature over a 
thermalization length, or length beyond which frequency 
redistribution affects the spectral flux profile 1111 . In this 
model the thermalization length is such that the final tem- 
perature is T = 4.0 x 10 4 K. The electron density n e is 
set to either 10 12 cirr 3 , 10 13 cirr 3 or 10 14 cirr 3 . Setting 
a constant n e implies that the velocity curve and/or the 
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FIGURE 1. Sobolev optical depth as a function of frequency 
at solar metallicity. 



stellar radius is different for each model, though neither 
as specified explicitly, as this study is concerned more 
with modeling a representative star, rather than a grid of 
specific stars. 



RESULTS 

Not surprisingly, the lack of metals, especially of iron 
group ions, causes a large drop in the line opacity. Figure 
[1] shows the Sobolev optical depth spectrum for a solar 
metallicity star with stellar parameters similar to those 
explored in this paper. Figure [2] shows the Sobolev op- 
tical depths for the Z = case in CRD. Clearly there is 
much less opacity to drive the wind. 

The mass loss rates and wind parameters are shown in 
table[T] Only models that produced physical results (e.g., 
Fred < 1) are shown. None of the CRD models were able 
to produce a wind. Thus for the given parameter ranges 
a CAK-type wind requires a finite thermalization length. 
As the wind density drops, M(t) rises, meaning when the 
continuum is weaker, the lines take on a larger fraction 
of the momentum driving. The models for which 10 < 
M{t) < 50 are most similar to Population I WR stars, and 
are probably the most plausible for Population III WR 
stars. Of these, the mass-loss rate is in the range 2 x 10~ 9 
M y r_1 < M < 3 x 10~ 8 M yr -1 . Compared to Vink 
and de Koter's low-Z mass-loss rate M w 10~ 7 M yr 1 , 
my rates are small, though the results are not directly 
comparable due to the difference in stellar mass and the 
WO-like abundances of models in this presentation. 



CONCLUSIONS 

A small parameter study of an analytic, modified CAK- 
type model of Population III WO-type stellar winds has 
been performed. The most plausible results were found 



for T = 0.5 and T = 0.7 and 10 13 cirr 3 < n e < 10 
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TABLE 1. Results for model calculations 



n e (cm J ) 


r 




t 




M(t) 


M (M Q yr x ) 


1 x 10 14 


0.3 


5 


3.6 x 10 


-5 


96.3 


7.3 x 10~ 9 


1 x 10 14 


0.3 


10 


4.1 x 10 


-5 


91.5 


8.6 x 10~ 9 


1 x 10 14 


0.5 


5 


1.0 x 10 


-4 


6.7 


1.6 x 10~ 8 


1 x 10 14 


0.5 


10 


1.1 x 10 


-4 


12.0 


1.9 x 10~ 8 


1 x 10 14 


0.7 


2 


5.0 x 10 


-4 


15.0 


1.9 x 10~ 8 


1 x 10 14 


0.7 


5 


7.1 x 10 


-4 


13.5 


2.5 x 10~ 8 


1 x 10 14 


0.7 


10 


8.9 x 10 


-4 


10.6 


2.8 x 10" 8 


1 x 10 13 


0.3 


5 


1.2 x 10" 


-5 


134 


8.6 x 10~ 10 


1 x 10 13 


0.3 


10 


1.7 x 10 


-5 


98.9 


1.0 x 10~ 9 


1 x 10 13 


0.5 


2 


5.8 x 10" 


-5 


49.7 


2.2 x 10~ 9 


1 x 10 13 


0.5 


5 


7.5 x 10" 


-5 


31.5 


3.5 x 10~ 9 


1 x 10 13 


0.5 


10 


8.6 x 10" 


-5 


26.4 


4.5 x 10~ 9 


1 x 10 13 


0.7 


2 


1.5 x 10" 


-4 


17.2 


6.5 x 10~ 9 


1 x 10 13 


0.7 


5 


1.7 x 10" 


-4 


15.8 


9.5 x 10~ 9 


1 x 10 13 


0.7 


10 


1.8 x 10" 


-4 


15.1 


1.2 x 10~ 8 


1 x 10 12 


0.3 


2 


2.8 x 10" 


-6 


377 


2.8 x 10~ 10 


1 x 10 12 


0.3 


5 


3.1 x 10" 


-6 


394 


3.9 x 10~ 10 


1 x 10 12 


0.3 


10 


3.3 x 10" 


-6 


326 


4.6 x 10~ 10 


1 x 10 12 


0.5 


2 


5.1 x 10" 


-6 


222 


7.2 x 10~ 10 


1 x 10 12 


0.5 


10 


5.8 x 10" 


-6 


225 


1.1 x 10~ 9 
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FIGURE 2. The same as figure[T] but at Z = and CRD. 



cm 
yr 



with mass-loss rates in the range 2 x 10 M fi 



-'<M<3x 1(T 8 M yr" 1 



A similar analysis could be done for WN stars, a phase 
through which all WR stars pass. WN stars will have an 
even lower mass-loss rate, as the only line driving would 
be provided by He and a relatively small amount of N. 
These results will be used in a study of mass-loss as a 
function of small but nonzero metallicity. 



ACKNOWLEDGMENTS 



Alamos National Security (LANS), LLC under contract 
No. DE-AC52-06NA25396. 



REFERENCES 

1. D. J. Hillier, and D. L. Miller, ApJ 496, 407-+ (1998). 

2. G. Grafener, and W.-R. Hamann, A&A 432, 633-645 

(2005) . 

3. J. I. Castor, D. C. Abbott, and R. I. Klein, ApJ 195, 
157-174 (1975). 

4. J. S. Vink, and A. de Koter, A&A 442, 587-596 (2005). 

5. V. V. Sobolev, Moving envelopes of stars, Cambridge: 
Harvard University Press, 1960. 

6. A. J. Onifer, and K. G. Gayley, ApJ 636, 1054-1066 

(2006) . 

7. K. G. Gayley, S. P. Owocki, and S. R. Cranmer, ApJ 442, 
296-310(1995). 

8. M. J. Barlow, and D. G. Hummer, "The WO Wolf-Rayet 
stars," in Wolf-Rayet Stars: Observations, Physics, 
Evolution, edited by C. W. H. de Loore, and A. J. Willis, 
1982, vol. 99 of IAU Symposium, pp. 387-392. 

9. R. L. Kingsburgh, M. J. Barlow, and P. J. Storey, A&A 
295, 75-100 (1995). 

10. N. R. Badnell, M. A. Bautista, K. Butler, F. Delahaye, 
C. Mendoza, P. Palmeri, C. J. Zeippen, and M. J. Seaton, 
MNRAS 360, 458-464 (2005). 

11. L. B. Lucy, and D. C. Abbott, ApJ 405, 738-746 (1993). 



The author would like to thank Alexander Heger for 
helpful discussions. This work was performed under the 
auspices of the U. S. Department of Energy by the Los 



